Combining Weak Lensing Tomography with Halo Clustering to Probe Dark Energy 
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Two methods of constraining the properties of dark energy are weak lensing tomography and 
cluster counting. Uncertainties in mass calibration of clusters can be reduced by using the properties 
of halo clustering (the clustering of clusters) . However, within a single survey, weak lensing and halo 
clustering probe the same density fluctuations. We explore the question of whether this information 
can be used twice - once in weak lensing and then again in halo clustering to calibrate cluster 
masses - or whether the combined dark energy constraints are weaker than the sum of the individual 
constraints. For a survey like the Dark Energy Survey (DES), we find that the cosmic shearing of 
source galaxies at high redshifts is indeed highly correlated with halo clustering at lower redshifts. 
Surprisingly, this correlation does not degrade cosmological constraints for a DES-like survey, and 
in fact, constraints are marginally improved since the correlations themselves act as additional 
observables. This considerably simplifies the analysis for a DES-like survey: when weak lensing 
and halo clustering are treated as independent experiments, the combined dark energy constraints 
(cluster counts included) are accurate if not slightly conservative. Our findings mirror those of 
Takada and Bridle, who investigated correlations between the cosmic shear and cluster counts. 
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I. INTRODUCTION 

Weak lensing tomography [E H, H, 0, H, Q and cluster 
counting 0, H, H, E3, [LH are complimentary methods of 
constraining dark energy parameters. Weak lensing to- 
mography (WLT) measures how the images of galaxies 
at various redshifts are sheared by large scale structure, 
while cluster counting (CC) simply counts the number 
of massive clusters as a function of redshift. Both tech- 
niques exploit dark energy's geometric effects and its ef- 
fects on large scale structure growth. The spatial cluster- 
ing of clusters, or halo clustering (HC), provides free ad- 
ditional information that can break degeneracies in clus- 
ter surveys [Til . [l2| . Although halo clustering is by far 
the weakest of the three probes, it becomes more valu- 
able as systematic uncertainties (e.g. cluster mass cali- 
bration) degrade dark energy constraints. Future surveys 
will combine these techniques. 

Since both the cosmic shear and the presence of clus- 
ters arise from fluctuations in the cosmic mass density 
field, we can expect their signals and error bars to be cor- 
related. Intuitively, if a survey or combination of surveys 
observes a patch of sky with more clusters than average, 
we should expect the patch to have an above average cos- 
mic shear signal. Hence, an analysis that treats lensing 
and clusters as independent observables will double-count 
information and may potentially overstate the bounds on 
cosmological parameters. On the other hand, the cross- 
correlation between lensing and clusters is an additional 
observable that can work to improve the parameter fit. 
Thus, combining lensing and clusters properly is impor- 
tant: treating them as independent could mean that we 
have artificially small error bars or that we are throwing 
away valuable information - both disturbing possibilities! 

Several authors have investigated WLT/CC correla- 



tions, finding that they arise primarily from the small 
scale power that individual halos contribute to the mass 
power spectrum. Approximating the halo number den- 
sity as a smooth field related linearly to the cosmic mass 
density, Fang and Haiman [13l ] find negligible covariance 
between cosmic shear and total cluster counts. Takada 



and Bridle [14| further take into account the structure 
and discrete nature of halos; they find a correlation be- 
tween cluster counts and the shear power spectrum on 
arcminute scales of about 10% (higher or lower depend- 
ing on the minimum halo mass used). 

We focus here on WLT/HC correlations: these should 
occur since shear maps and halos trace the same field of 
projected matter fluctuations (on scales larger than the 
average halo spacing). In other words, a local overdensity 
of mass tends to create a shear peak and an abundance 
of halos along the same line of sight. This correlation 
turns out to be significantly larger than the WLT/CC 
correlation, and it warrants a closer look at how halo 
clustering and weak lensing tomography should be com- 
bined. The WLT/HC correlation is most important in 
the regime where halo clustering is adding significant con- 
straints to dark energy parameters; we will be working 
in that regime. 

In this paper, we calculate the WLT/HC cross- 
correlation and determine its effect on a survey's ability 
to measure dark energy parameters. We restrict ourselves 
to large angular scales so that halo clustering is related 
to the mass overdensity by a linear bias and we can ig- 
nore the effects of individual halos on the shear power 
spectra. We focus on projections for the Dark Energy 
Survey (DES) and find that including the full covariance 
does not degrade the projected constraints; in fact these 
tighten at the percent level. 
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II. WEAK LENSING AND HALO CLUSTERING 

Divide a survey region into redshift bins with Zi denot- 
ing the minimum redshift of the ith bin. In this paper, 
we consider 10 bins of equal width up to z=l. Let x( z ) 
denote the comoving distance to an object at redshift 
z so that Xi = x( z i) (we will often use x an d z inter- 
changeably). The weak lensing convergence is a good 
computational proxy for shear; in a given redshift bin, 
the convergence is a function of angular position on the 
sky and can be expressed as a weighted integral of the 
density along the line of sight: 



where the weighting function is 
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for x < Xi+i an d otherwise [15J. Here vf al is the ex- 
pected galaxy density in bin i which we parametrize as 



(1) 



dv^ 1 
dXs 
(2) 



dv 



gal 



dz 



-z I Z* 



(3) 



The distribution has z mo de = 2-z, and z mc dian = 2.67z* ; 
for DES we use z* = 0.254 and ^f al = 15/arcmin 2 . 

We approximate the local overdensity of halos in the 
ith redshift bin as another line-of-sight integral, 



hi(0) 



_ 6m(8) 



dxW (Hai oMx)S[xO,x}, (4) 



where fii is the average halo density in bin i. We are 
counting all halos more massive than a threshold M m j n = 
10 14 - 2 M Q //i, and we've assumed that the overdensity of 
halos of a given mass is proportional to S (valid on large 
scales). The halo weighting function is 



W, 



(Halo,») 



X 2 b M (x) Xi<X< Xi+i (5) 
otherwise 



where the mass-averaged halo bias b M {x) depends on 
both redshift and our M m - ln : 



W( , lM ml JMb(M,z)^(z) 



(6) 



We compute b M using the halo mass function dn/dM(z) 
of Jenkins et al. flq| and the halo bias model b(M, z) of 
Sheth and Tormen \v\. 

Let ui be one of our observables with the subscript 
denoting the data type and redshift bin, e.g. uj — for 
/ = (Lens, i). Then 



Ul (0) = J d x w I (x)S[xe, x }- 



(7) 
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FIG. 1: Sample correlation coefficients; for a given lensing 
bin, the most strongly correlated halo bin is shown. All bins 
have width Az = 0.1, and the average redshifts of the bins 
(Lens, Halo) are as follows: solid blue (0.95, 0.25); dot-dashed 
green (0.65, 0.15); dashed purple (0.35, 0.05). 



Using the flat-sky approximation, decompose the observ- 
ables uj into 2D Fourier modes. Then we define 



( Ul (\)uj(l')) = (2nf6 2 (l-l')C I j(l) 
where under the Limber approximation, 
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The matter power spectrum Pu(k) depends on the data 
types. We use a linear spectrum to compute halo-halo 
and halo-lens terms and a nonlinear spectrum to com- 
pute lens-lens terms 1 . Our matter power spectra are 
computed using the fitting formula of Eisenstein and Hu 
[l8| and the nonlinear scaling of Smith et al. [l9j . Note 
that the weighting functions Wi imply that 



(hihj) = for i ^ j 
(nihj) = for i < j. 



(10) 
(11) 



These follow from the Limber approximation and should 
be valid on scales I > 10. 

To assess the degree to which lensing and halo cluster- 
ing are correlated, it is useful to consider the correlation 
coefficient 



R U EE C U {l)/^C U {l)Cjj{l). 



(12) 



If lensing were completely uncorrelated with halo clus- 
tering, the correlation coefficients i?(Lons.i)(Haioj) would 
vanish; if they measured identical quantities, the coeffi- 
cient would be unity. Figure [1] shows that on large scales, 



1 In the formalism we've chosen, the proper treatment of the halo- 
lens term is ambiguous, but our results are insensitive to our 
matter power spectrum choice. 
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the lensing of high redshift background galaxies is highly 
correlated with clustering at intermediate redshifts. It 
is therefore possible that tighter constraints anticipated 
from including halo clustering are illusory: perhaps this 
information has already been used once - in weak lens- 
ing tomography - and can't be double-counted. In that 
case, it would not be appropriate to blindly combine the 
constraints from (CC+HC) with those from WLT. 



as Gaussian): 

where <5/j is a Kronecker delta and 

7^/i/f ;al for I = (Lens, i) 



N I = 



1/rii for 7 = (Halo, i) 



(15) 



(16) 



III. IMPACT ON DARK ENERGY 
PROJECTIONS 

To assess the impact of the correlations depicted in 
Fig. [TJ we compute the Fisher matrix for weak lensing 
tomography and cluster counting supplmented by halo 
clustering. We choose a standard set of cosmological pa- 
rameters in a flat universe: 



Parameter 




n b h 2 


Side U)0 


w a n 


*c 2 


Fiducial 


.14 


.024 


.73 -1 


1 


2.57 x 1(T 9 


Prior 


1% 


1% 




- 1% 


1% 



where h is the Hubble constant in units of 100 
km/sec/Mpc; fide is the ratio of the dark energy den- 
sity to the critical density (and similarly for the total 
matter and baryon densities); n is the spectral index of 
primordial scalar fluctuations, and S( is their amplitude 2 
at k = 0.002 Mpc -1 ; the dark energy equation of state is 
taken to be w — wo + w a z/(l + z). We impose Gaussian 
priors expected from future CMB surveys. Additionally, 
we include (with no priors) the cluster mass calibration 
parameters A and tia defined by Lima and Hu [l2j ■ 

The Fisher matrix for our parameter set is, in the stan- 
dard approximation, a sum of matrices for our three ob- 
servables: 
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(13) 



The full cluster Fisher matrix breaks nicely into a CC 
and HC term, as in the above equation [12|]. However, 
we know from the results of §11 (e.g. Fig. [1]) and pre- 
vious work ([H [l|) that lensing and clusters are not 
completely independent observables, and therefore, that 
the above approximation is flawed. To account for the 
correlation between halo clustering and cosmic shear, we 
generalize F Total to 



TiTotal 
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CC 



,-iHC+WLT 



Ti WLT 



(14) 



where F WLrr contains lensing information on scales 
smaller than can be probed with halo clustering. 

The observed spectra C ohs (l) are the sum of the sig- 
nals given by Eq. © and Poisson noise (approximated 



Here, 7* is the RMS intrinsic shear of the source galaxies; 
for a DES-like survey, 7* w 0.25 [2fJ. Let Q be a matrix 
with elements Cjj s (Z). Then for a set of parameters p a , 



pHC+WLT 
c a(5 



/sky 
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I 
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dQ c _ t 8Q c _ x 
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(17) 

where / s k y is the fraction of sky covered by the survey, 
and the trace sums over all redshift bins and data types 3 . 
Similarly, we can write i^ WLT this way, summing only 



over the lensing part of Q. 
sum over 



HC+WLT 



is a 



Note that F 
while _F WLT is summed 



from L 



T to L 



We choose L 



=20 and L 



=200 



so that the flat-sky approximation and our linear halo 
bias assumption are valid. As prescribed by Rudd et al. 
[2TI ] . we modestly set / max =1000 since baryonic effects on 
smaller scales are not completely understood. 

If halo clustering were a completely independent ob- 
servable from cosmic shear, then the matrices Q would 
be block diagonal with no terms in the (HC,WLT) re- 
gion. In that case, Eq. (fTl)) reduces to Eq. (fT3"|) . We 
compute i^ Total for both cases, and since we are not in- 
vestigating CC correlations, our cluster counting term is 
always given by 



P cc 



x - drrii drrii 1 
4-* dp a dpp mi 



(18) 



where m, is the total number of clusters in the zth redshift 
bin with M > M m i n . Note we've assumed that Poisson 
noise is the dominant source of scatter for cluster count- 
ing. 

After computing the Fisher matrix, we use it to fore- 
cast la errors for the dark energy parameters after 
marginalizing over all other parameters. We also com- 
pute the dark energy figure of merit (FoM), defined as 
the inverse area of the 95.4% error ellipse for (wo,w a ) 
[22I ]. It turns out that the HC/WLT correlations have 
very little impact on the final parameter determination 
and even bring a slight improvement for DES. Numbers 
quantifying this small impact are given in Fig. [2] and Ta- 
ble |U We remind the reader that Fig. [5] varies two of 



3 We check that our numerical derivatives converge with decreasing 
step size. Ultimately, we use two-sided derivatives with steps that 
2 In practice, we vary ln8 2 ( . The prior is defined at k = are 4% of the paramete rs or 0.04, 0.01 and 0.01 for w a , A and 

0.05 Mpc -1 . respectively. 
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FIG. 2: Fractional change in the dark energy figure of merit 
once WLT+HC covariance is included, shown as a function of 
sky coverage and galaxy density. Positive values imply that 
neglecting the covariance results in an underestimation of the 
FoM (conservative parameter constraints). 



Parameter 


No HC 


WLT+HC 


WLT+HC 


Difference 






Correlated 


Independent 




a(wo) 


0.601 


0.461 


0.470 


-2% 


<?{w a ) 


2.15 


1.52 


1.56 


-2% 


cr(fidc) 


0.0387 


0.0362 


0.0377 


-4% 


FoM 


0.485 


1.15 


1.01 


+ 14% 



TABLE I: la error forecasts and dark energy figure of merit 
for a DES-like survey (/ s k y =0.12, j/f al =15/arcmin -2 ). Con- 
straints in the 2nd column combine CC+WLT only. Columns 
3 and 4 show constraints from CC+WLT+HC, treating 
WLT+HC as correlated or as independent observables. The 
last column compares columns 3 and 4. 



several possible survey parameters: M m - m , z* and 7, are 
held fixed. 

Our results show that cosmological constraints can be 
marginally improved by including the full HC/WLT co- 



variance since the HC/WLT correlation itself acts as 
an additional observable. Whatever degradation occurs 
from acknowledging that lcnsing and clusters are corre- 
lated may be offset by our ability to predict and measure 
the correlation. We note that correlations do not nec- 
essarily imply error degradations in the first place: this 
would be the case if our set of observables were correlated 
in such a way that does not mimic their dependence on 
dark energy parameters. Mirroring the results of Takada 
and Bridle [14| , we conclude that a DES-like survey that 
treats halo clustering and weak lensing tomography as in- 
dependent experiments will have accurate if not slightly 
conservative parameter constraints. 



IV. CONCLUSIONS 

The distribution of halos along a given line of sight 
probes the mass density along that line of sight in ways 
very similar to weak lensing. We have quantified this 
similarlity with the cross-correlation coefficients shown 
in Fig.Q] However, for a survey like DES, projected con- 
straints on dark energy parameters are not degraded by 
this cross-correlation (using the same information twice), 
and it is acceptable to simply add the constraints from 
weak lensing to those from clusters. Doing so actually 
results in slightly conservative parameter constraints. 
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